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Abstract. - The study of dynamo action in astrophysical objects classically involves two 
timescales: the slow diffusive one and the fast advective one. We investigate the possibility 
of field amplification on an intermediate timescale associated with time dependent modulations of 
the flow. We consider a simple steady configuration for which dynamo action is not realised. We 
study the effect of time dependent perturbations of the flow. We show that some vanishing low 
frequency perturbations can yield exponential growth of the magnetic fleld on the typical time 
scale of oscillation. The dynamo mechanism relies here on a parametric instability associated with 
transient ampliflcation by shear flows. Consequences on natural dynamos are discussed. 
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Introduction. - Most astrophysical objects exhibit mag- 
netic fields. This is the case of planets, stars, as well 
as galaxies. The dynamo instabifity, characterised by 
an exponential growth of the magnetic field in a mag- 
netohydrodynamic flow, is expected to account for these 
fields. Among dynamo flows, a vast majority provides 
finite growth rate on the diffusion time only (known as 
''slow dynamos" ) . They are therefore not suitable for large 
astrophysical objects (stars and galaxies), for which the 
associated magnetic Reynolds number (measuring the rel- 
ative strength of the induction to the diffusive term) is 
huge. This is the primary motivation for the investiga- 
tion of fast dynamos, characterised by finite growth rate 
on the advective time scale rather than the diffusive time 
scale (the later being sometimes larger than the age of the 
universe) . 

While significant progress has been achieved in the case 
of moderate values of the magnetic Reynolds number [1,2], 
the fast dynamo limit has been tackled yet only with a few 
analytic flows [3,4], still remote from realistic configura- 
tions. Criteria to assess dynamo action (other than direct 
integration) are still missing. Most simple flows (such as 
uniform shear) fail to be dynamos and only yield transient 
amplification of the field, followed by ohmic decay. 

We discuss the possibility of a third class of dynamo, for 
which the field grows on an intermediate timescale associ- 



ated to flow modulations (for example due to a traveling 
wave). The efficiency of some time dependent flows in 
the dynamo process has been recognised in various stud- 
ies [4-8]. These studies focused on leading order time 
dependence of the flow (either analytical, or in the form 
of an heteroclinic cycle), whereas we want to study the 
effects of vanishing perturbations on steady flows. We 
stress that the impact of small fluctuations on dynamo 
flows has already been addressed in other contexts. Gog 
et al. [8] considered the action of a (random) noise on the 
dynamo properties of time dependent flows provided by 
heteroclinic cycles. Petrelis & Fauve [9] and Peyrot et 
al. [10] investigated the effects of fluctuations on the dy- 
namo threshold (respectively for the G.O. Roberts and the 
Ponomarenko flows). 

The magnetic field evolution in a conducting fluid is 
governed by the induction equation 



^tB = V X (u X B) -f ryAB. 



(1) 



The rhs operator is not self-adjoint [11-14] and can there- 
fore provide algebraic growth of an initial perturbation 
[15]. This transient magnetic amplification is typically ob- 
served in a uniform shear (as examplified by the Jl-effect) . 
We will see how this effect can be used to yield dynamo in- 
stability under appropriate time dependent disturbances. 
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Modelling. - We focus on a very simple setting, re- 
taining only the essential physical characteristics of the 
above system (the transient amplification). Let us con- 
sider a two-dimensional flow u = {ux{y,t), 0, 0)* and fm-- 
ther assume an initial magnetic field of the form Bq = 
{Bx cos(z), By cos(z), 0)* independent of the position x. 
The field depends only on the space coordinate z, which 
provides a typical length scale L. 

We assume that the flow consists of a uniform and con- 
stant shear dyUx{y,t) = S . The evolution of Bx and By 
in time is then governed by 



dt 
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The time variation of the field strength is measured by 
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Bl + Bl dz 



{Bx{Q)+StBy{Q)Y + By{Qf 



(3) 



which clearly highlights the transient algebraic amplifica- 
tion for an initial seed field along By. 

We will now discuss the effect of a time dependent per- 
turbation of the flow. Two typical time scales can be iden- 
tified in the original problem (2). The advective time scale 
T„ = L/U = and the diffusive time scale t,, = i^/?7. 
As we focus here on the limit of large magnetic Reynolds 
numbers (relevant to astrophysical dynamos), we shall 
preferably use the advective time scale. The system can 
be rescalcd accordingly using e ~ rj/UL (the inverse of the 
magnetic Reynolds number), and 



-e 1 
-e 



(4) 



We will now consider small amplitude, time-periodic, 
perturbations of the flow u. These could for example 
represent the effect of a traveling wave. For simplicity, 
we assume that the perturbation modifies the direction 
of the flow but not its amplitude. The pulsation is de- 
scribed by a sinusoidal wave of amplitude a and frequency 
(3. The new velocity field can therefore be expressed as 
^ncw ^ 'jz{0{t))u , where 71(9) is the 2x2 rotation ma- 
trix of angle 9, and 9{t) = a sin(/3t). The evolution of the 
magnetic field B is then governed by the linear operator 



c{t)^n{9{t)) £„ n{9{t)y 



(5) 



Note that we have just introduced a third time scale in 
the problem, Tp ~ L/{f3U). 

We will investigate the field evolution and the possibil- 
ity of dynamo instability, yet the system was made very 
simple: the field is independent of x and y, and the fiow is 
planar [u^ = 0). It is then natural to ponder the feasibil- 
ity of dynamo action. In fact, the simplicity of the model 
will also prevent the application of the "two-dimensional 
field" and "planar velocity" anti-dynamo theorems [16,17]. 



As the field is independent of x and y, an expression in 
the form of a streamfunction would necessarily yield un- 
bounded values: the field considered here does not decay 
to infinity. The above anti-dynamo theorems therefore do 
not apply to this setup. 

Rewriting this system in the rotating frame relative to 
which the fiow is steady yields 



dt 



Bx 
B,, 



-e 1 + 9{t)\ (Bx 
-m -e ) \By 



(6) 



This problem can be tackled through the change of un- 
known B' = e^* B , and a rescaling in time to the unit 
of time, t' = (3t which leads to 



1//3 + a cos t' 
-a cos t' 



(7) 



Let us stress that, for simplicity, we have introduced this 
system in a very restricted setup, however such formula- 
tions can be extended to local behaviour in more complex 
fiows, for instance near stagnation points of pulsed Bel- 
trami waves [7] . 

Floquet analysis. - The stability of this system can be 
investigated depending on the values of a, (3. We will in 
particular focus on the limit of small a, i.e. the limit of 
infinitesimal perturbations of the fiow. This problem fits 
the framework of the Floquet analysis, and the stability 
is determined by the logarithm of the resolvent matrix 
eigenvalues. Numerical computation of these eigenvalues 
are reported Figure 1. For a perturbation of a definite 
temporal form (3, finite growth rate can only be achieved 
for finite amplitude modulations: there exists a thresh- 
old in the amplitude a below which a' vanishes (and the 
unfiltered variable B decays). This reflects the impossi- 
bility of fast dynamo action (i.e. growth on the advective 
time scale) with the chosen flow. However, in the limit of 
timescales separation, that is 



TU <. Tp <^ Tn , 



(8) 



which yields 0, exponential growth can occur on the 
Tp time scale, as revealed by the finite growth rate a' at 
any value of a for small enough (3. Tiny disturbances 
can yield exponentially growing solutions on a typical time 
scale provided by the perturbation (i.e. Tp). This yields 
the notion of time scales separation. Remarkably, the 
growthrate does not always vanish in the limit of small 
a (i.e. infinitesimal perturbation). 

Figure 1 reveals a scaling of the form a on (3, which 
points to the use of 7 = a//?. We investigate on Figure 2a 
the growth rate of (7) as a function of (a, 7). The system 
has an intricate behavior, with large regions of instability, 
even in the limit a ^ , at finite 7, separated by narrower 
neutrally stable bands. Finite growthrate on the T/3 time 
scale for large 7 (small /3 at fixed a) is even clearer than on 
Figure 1. This behavior is reminiscent of the parametric 
instability that develops in forced mechanical systems. A 
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Fig. 1: Growth rate a' (in the time scale) as a function of 
a and (3. 

typical example is provided by the Mathieu equation [18] 
(governing a vertically excited pendulum), 

6i+ (5 + 7cos<) 61 0. (9) 

The classical stability diagram for this system is repre- 
sented on Figure 2b. Note that system (7) can be rewritten 
in the scalar form 

Q/-^ sin ij 

B'x^ 5 B' +cost(7 + a2 cosi)B' = 0. (10) 

7 + a cos t 

It is worth stressing that this system does not reduce in an 
appropriate limit to the unforced pendulum, i.e. 7 = in 
(9). An equivalence between both problems can, however, 
be estabhshed in the limit a, 5 (which corresponds to 
vanishing gravity for the pendulum). This is highlighted 
on Figure 2c which presents cross sections of Figures 2a 
and 2b for a = and 5 = respectively. 

The growth rate (T'(a,7) seems unbounded as 7 —> cx), 
but this is only an effect of the chosen time scale , which 
becomes increasingly long compared to t„ in this limit. 
Translating to the advective time scale Tu, the growth 
rate is given by a{a,^) = /3a'{a,j) = \a\a' {a,j)/\j\. 
The quantity cr" = cr' (a, 7)/ I7I is represented on Fig- 
ure 3, which clearly shows that cr" is uniformly bounded 
for all (0,7) and reaches positive values. This validates 
the scaling a ^ (3, the optimum growth rate is achieved 
for 7 « 0.786, and behaves like a « 0.454 a for small a. 

One should keep in mind that ohmic decay has been fil- 
tered out from (6) to (7) : exponential growth therefore re- 
quires a > £ and thus a > e. This corresponds to ensuring 
that a modified magnetic Reynolds number constructed as 
Tri/rp remains larger than unity. 

Interpretation. - It should be stressed that the insta- 
bilities reported here are only made possible by the non 
self-adjoint nature of the C operator, i.e. of the induction 
equation. If the system was governed by a negative self- 
adjoint operator, it could be readily shown that a time 
dependence of the form (5) would imply negative values 
of 

i dt|B|2 = (/:7^-l(0(^)) b | Ti-\e(t)) b) < o . (ii) 
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Fig. 2: (Color online) a. Growth rate cr' represented in the 
(7, a)-plane. b. Growth rate for the Mathieu equation in the 
(7, 5)-plane. c. Both solutions are identical when a = 5 = 0. 

The role of transient growth in this parametric insta- 
bility can be further analysed by taking advantage of the 
fact /? ^ 0. The system can then be investigated under 
the quasi-static approximation of constant = ao, and 
locally replaced near each time by 

C{t)^n{aat) Cun{aQt)-\ (12) 

where = a(3 cos(/3to). This system corresponds to a 
steady rotation of the shearing flow at angular velocity ao 
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Fig. 3: (Color online) Growth rate cr" = a' [a,^) in the 
(7, a)-plane. Maxima are achieved at finite values of 7. 

(such situation can occur in convcctivcly driven dynamos, 
e.g. [20]). It is completely integrable, and admits expo- 
nentially growing solutions for —1 < ag < , 

o- = -£ + V-ao (1 + Q'o) • (13) 

This formula is reminiscent of the classical Rayleigh- 
Pedley dispersion relation for rotating shear flows, in 
which ao is replaced by 2 [19]. This similarity is due to 
the well known formal analogy of the induction and vor- 
ticity equations. The governing pertubation equations do 
however differ. 

This exponential growth can be directly related to 
the transient field amplification in the initial system (2). 
While an initial magnetic field is amplified by the shear as 
described by (3), its direction is altered. This results in 
a decreases of the phase with time. Let us consider this 
mechanism over a small increment of time 5t. The initial 
field is amplified 

\m + 5t)\\=e^''\\B{t)\\. (14) 
and rotated by a small angle 

B(t + 5i)/||B(t + (5t)|| =7^(S,5^)B(^)/||B(^)||. (15) 

If we now consider a system of the form (12), ao can 
be set to a so as to cancel the effect of rotation, which 
leaves us with 'R{t-\-5t) = e°'''* B(t). Incrementing in time 
then yields an exponential increase of the amplitude with 
growth rate a. Initially, the field which is the most am- 
plified by (2) is (1, Vf I\f2, and corresponds to tr = i — e 
and a = — ^. This is precisely the maximum growth rate 
in (13), and establishes the direct connection between the 
change in the field orientation during transient amplifica- 
tion and the exponential growth obtained with (12). The 
rotation rate a. provides a typical time scale oT^ for sys- 
tems exhibiting algebraic growth. This time scale will pro- 
vide the physical interpretation of the instability process 
investigated here. 
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Fig. 4: Growth rate in the advective time scale versus a/3. 
Dashed lines correspond to the result of a Floquet analysis on 
system (7) at fixed values of a, from right to left respectively 
a = 0.25,0.1875,0.125,0.0625. The solid black curve corre- 
sponds to the estimate (16). 

If we return to (6) and use the approximation (12), the 
rotation rate is provided by a/3cos(/3to)- For small enough 
a/?, formula (13) can be used over half a period, corre- 
sponding to cos(/3fo) < 0. This expression being related to 
the fastest growing mode, the phase of the solution will act 
to reduce the realized growth rate. During the other half 
of the period (cos(/?to) > 0) , algebraic growth is achieved. 
A good approximation of the resulting growth rate can be 
obtained by simply considering the growth rate provided 
by (13) and ao = a/3 on the first half period and a zero 
growth rate for the reminder. As a/3 ^ 1 this yields at 
leading order 

a = \^fafi. (16) 

Direct comparison of this expression (still on the advective 
time scale) with the Floquet integration is provided on 
Figure 4. 

Expression (13) also sheds light on the requirement e —* 
for instability in the limit of infinitesimal perturbations. 
Indeed a/3 > is needed to allow positive values of tr. 

The exponential instabilities investigated here can 
therefore be understood both in terms of parametric insta- 
bilities and in terms of perturbed algebraic amplification. 
It is interesting to note that the steady operator L is not 
an oscillator and has no characteristic frequency. Since it 
is associated with transient amplification, a typical time 
scale can, however, be defined as introduced above. 
The evolution of the unperturbed system provides, as for 
the pendulum (Mathieu equation), the natural time scale 
for parametric amplification. 

Conclusions. - In the limit of large magnetic Reynolds 
number Rm = e~^, a distinction is usually made between 
dynamo action on the "slow" diffusive time scale and 
the "fast" advective time scale r^. We stress here the po- 
tential importance of a third time scale corresponding to 
flow modulations (here r^), on which exponential growth 
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can be achieved. Natural dynamos often exhibit large scale 
shear (usually in the form of a differential rotation associ- 
ated with the so-called "fi-effect" ) . Such large scale shear 
are unable to produce dynamo action on their own. Al- 
though the model we present is extremely simplified (the 
shear is uniform and the field only depends on one spatial 
coordinate), we believe the conclusion that interaction of 
shear and rapid fluctuations can result in dynamo action is 
an important property. Although we restricted our atten- 
tion to a simplified setup, interesting conclusions can be 
drawn concerning dynamo action. In the limit of large i?m, 
the flow needs to be measured with a precision Rm~^ for 
its dynamo property to be determined. Tiny fluctuations 
(as small as Rm~^) can drastically change the dynamo 
properties of the flow. In astrophysical objects RmT^ is 
often as small as 10^^*, so tiny fluctuations (but larger 
than RrrT^) can modify a non-dynamo flow to provide ex- 
ponential growth. Further studies are obviously needed 
to assess under which conditions this amplifications is re- 
alised in fully three dimensional configurations. 

The possibility to drive a dynamo through periodic 
modifications of the flow has already been pointed out, 
but such an approach was so far restricted to finite am- 
plitude modulations. We find that even small amplitude 
oscillations can act to drive parametric instabilities, while 
the unperturbed operator is stable. More generally, many 
physical problems are governed by non self-adjoint oper- 
ators. The system we investigated being very simplified, 
it could be interpreted in a much more general framework 
than that of the induction equation (for instance in hy- 
drodynamics, see [15,21]). 
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